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ABSTRACT 

A  circular  cylinder  immersed  in  an  infinite  fluid  and  oscillated 
perpendicular  to  its  axis  is  acted  upon  by  a  fluid  dynamic  force.   This 
force  is  generally  considered  to  be  composed  of  two  components,  a  drag 
and  an  added  mass  component.   The  relative  importance  of  these  are 
functions  of  the  frequency  and  displacement  of  the  cylinder. 

It  was  the  purpose  of  this  study  to  determine  experimentally  these 
two  components  as  a  function  of  the  ratio  of  the  displacement  amplitude 
to  the  cylinder  diameter,  and  the  Reynolds  number  based  on  the  maximum 
velocity. 
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I.   INTRODUCTION 

With  the  increased  deployment  of  submerged  ocean  structures  the 
need  to  predict  the  forces  exerted  on  bodies  immersed  in  an  oscillating 
flow  is  of  practical  importance.   Studies  conducted  by  Reid  and 
Bretschneider  [1]  on  submerged  structures  stress  the  necessity  of 
having  these  forces  forecast,  given  a  stated  set  of  wave  parameters. 

A  number  of  studies  have  been  carried  out  to  determine  the  values 
of  the  added  mass  and  drag  coefficient  for  a  pile  in  ocean  waves.  Many 
of  these  are  represented  in  a  correlation  of  drag  coefficient  versus 
Reynolds  number  presented  by  Wiegel  [2].   However,  his  plot  shows  an 
order  of  magnitude  scatter  in  the  values. 

Few  basic  studies  of  forces  exerted  on  bodies  in  oscillating  flow 
have  been  conducted.   Of  these,  however,  the  most  well  known  was 
carried  out  by  Keulegan  and  Carpenter  [3].   They  conducted  tests  with 
a  stationary  cylinder  placed  at  the  node  of  a  standing  wave  and  obtained 
a  correlation  of  Cm  (mass  coefficient)  and  Cd  (drag  coefficient)  versus 
Um  T/D,  where  Um  denotes  the  maximum  velocity,  T  the  period,  and  D  the 
cylinder  diameter.   It  was  suggested  by  Keulegan  and  Carpenter  that  the 
parameter,  Um  T/D,  or  equivalently ,  the  relative  fluid  displacement^ A/D, 
was  the  parameter  of  primary  importance;  the  effect  of  Reynolds  number 
was  assumed  to  be  of  little  influence  on  the  values  of  the  drag  and  added 
mass  coefficients.   It  is  noted,  however,  that  the  test  method  used  by 
Keulegan  and  Carpenter  allowed  no  control  over  the  Reynolds  number  so  that 
they  had  no  convenient  method  of  testing  this  influence. 


Heinzer  [4]  conducted  visual  studies  of  an  oscillating  circular 
cylinder  and  concluded  that  the  wake  characteristics  were  primarily 
dependent  upon  amplitude  of  oscillation  vice  Reynolds  number.         i 
Sarpkaya  and  Garrison  [5]  investigated  unidirectional  flow  past  a 
cylinder  with  constant  acceleration.   Their  work  lead  to  the  fact  that 
drag  and  inertial  forces  are  representable  as  functions  of  the  relative 

displacement  of  the  fluid.   These  studies  were  conducted  in  the  range 

4 
of  Reynolds  numbers  near   10   and  no  attempt  was  made  to  investigate 

the  effect  of  the  Reynolds  number. 

In  the  present  study  various  circular  cylinders  were  oscillated  with 

simple  harmonic  motion  in  water  at  rest.   Both  the  amplitude  and  the 

frequency  of  oscillation  were  controllable  so  that  both  the  effect  of 

Reynolds  number  and  relative  displacement  of  the  fluid  could  be  studied. 

However,  time  limitations  restricted  the  experimental  program  to  a 

series  of  runs  carried  out  at  one  fixed  Reynolds  number  with  increasing 

relative  displacement.   Also,  at  a  fixed  value  of  relative  displacement 

the  effect  of  Reynolds  number  was  determined. 


II.   METHOD  OF  APPROACH 

This  study  utilized  test  cylinders  of  4  inches,  2.75  inches  and 
1.5  inches  diameter,  16.5  inches  long  mounted  on  a  dynamometer.   The 
dynamometer  was  basically  composed  of  a  cantilever  beam  affixed  with 
strain  gauges.   The  strain  gauge  bridge  error  was  read  out  on  a  strip 
chart  recorder  calibrated  to  read  directly  in  pounds  force. 

The  force  traces  were  read  visually  at  twenty  intervals  over  a 
complete  cycle  and  this  numerical  force  data  used  as  an  input  to  a  data 
reduction  program  (Appendix  A  and  B)  based  on  Morison's  equation. 

Morison  et.  al.  [6]  consolidated  the  three  part  equation  of  McNown 
and  Wolf  [7]. 

F  =  Ao  p   ^- ^y   +  5  Px  ds  +  J  Cd  D  p   U|U|  (1) 

by  combining  the  first  two  terms  to  yield 

F  =  Cm  p  Ao  ^  +  Y  Cd  D  p   U|U|  (2) 

where     j-  Force  per  Unit  Length 

Ao  Circular  Cross  Section  Area  of  Cylinder 

D  Body  Dimens  ion  Normal  to  Flow 

p   Density  of  Fluid 

t   Time 

Cd  Coefficient  of  Drag 

Cm  Added  Mass  (Inertia)  Coefficient 

U  Velocity  at  Time  t 


Equation  (2)  is  assumed  valid  and  is  used  as  the  basis  of  the  data 
reduction. 

In  the  present  instance  the  velocity  of  the  cylinder  is  represented 
by: 

U  =  -Urn  cos  a  t  (3) 

where  Um  denotes  the  maximum  velocity,  T  the  period,  and  a  =  27t/T. 
The  total  force  acting  on  the  cylinder  per  unit  length  is  in  general 
given  by 

F  =  f(t,  T,  Um,  D,  P,  V)  (4) 

Grouping  the  variables  on  the  basis  of  dimentional  reasoning  and 
introducing  a  =  2iTt/T  gives 

F     ^,        UmT  UmD.  ,^. 

=  f(c,  -;r-,  — )  (5) 


Um  D 

,      UmP   .    „    -ij     1-     J  UmT    ,  .  ,       ,    , 
where  -— -  xs  a  Reynolds  number  and  "pr"*  which  can  also  be 

2TrA 
expressed  as  ~~^i     will  be  termed  the  relative  displacement. 

Because  of  flow  symmetry  and  the  periodic  nature  of  the  force 

F(a)  =  -F(a  +  tt)  (6) 

it  is  possible  to  express  the  force  coefficient  in  a  Fourier  series, 
F 


„  =  A^  sin  a  +  A»  sin  3a  +  A^  sin  5a  +  . . . 
Um  D  ^  ^ 


+  B^  cos  a  +  B»  cos  3a  +  B-  cos  5a  +  . . . 


(7) 
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where  the  coefficients  A^,   B   are  independent  of  a  and  at  most 
functions  of  Reynolds  number  and  relative  displacement.   Fourier  analysis 
may  be  used  to  determine  the  coefficients  as:  \ 


■2-n 
A     ■•-  \    r  sin  JNg  ,  .-v 

"     0     pUm  D 


1  \    F  sxn  Na 


and 


.2Tr 


0     pUm  D 


Once  obtained,  the  dependence  of  these  coefficients  on  Reynolds  number 
and  relative  displacement  may  be  established  provided  the  data  are 
sufficient. 

The  general  formulation,  equation  (7),  may  be  reconciled  with 
Morison's  equation  (2).   Introducing  U  from  equation  (3)  into 
equation  (2) 


Cm  •  TT-   sin  a  -  — r   cos  a   cos  a  (10) 


,.  2_   4      Urn  ^^"  "    2 
pUm  D 


By  the  rule  of  Fourier 

■  2tt 
cos  a  I  cos  a  =  ^   0 


■^  cos  a   cos  a  cos  Nada 


N=0  C^""  2 

J         cos  Nada 

0 


=  a^  +  a^  cos  a  +  a«  cos  2a  +  a_  cos  3a  +  ... 


where 

N+1 


8 


a     =  0,   a  ^^  =  (-1)   2  _o (u) 

even        odd   ^   '     -,,„2  ,.  , 

N(N  -4)tt 
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The  first  three  non-zero  coefficients  would  then  be: 


a  =  8/37r,    a  =  8/157r,    a^  =  S/lOSu  (12) 


Introducing  equation  (12)  into  equation  (7)  with 


h   '  h'^i 


h  ==3-V^iPi, 


t 
'5  "  "5       "5' 


B,  =  B^  -  a^/a^  (bJ 


(13) 


yields 


F 

=  A-  sin  a  +  A_  sin  3a  +  A^  sin  5a  +  . . . 


..  2^   "1  ^^"  ^       "3  "^"  ^^    "5 
pUm  D 

f 

-»■  B-   I  cos  a  I  cos  a  +  B„   cos  3a  +  B_   cos  5a  ■♦-..- 


(lA) 


Equations  (14)  and  (7)  may  be  compared.   Writing 

f  Cm  •  ^=  A,  +  A,  4^^-^+  A,  42-^+  ...  (15) 

4      Um    1    3  sin  a     5  sin  a 


and 


_,      I    B„   cos  3a  Bj.   cos  a 

^  =  _B   -  ^ , -  -p-2 +  . .  .      (16) 

2      1     I  cos  a  I  cos  a      |  cos  a  |  cos  a 

I      t 

Thus,  if  the  coefficients  A_,   A^.   and  B_  ,   B_   vanish,  the 

values  of  Cm  and  Cd  remain  constant  for  all  phases  of  the  cylinder 
motion 

n  2  UmT  .  .T,v 

Cm  =  -y  -3-  A^  (17) 

IT 
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Substituting  A   from  equation  (8)  yields 


•27T 


2   UmT  C^"   F  sin  a  da  .,.. 

Cm  =  —  -^  J    :i (18) 


TT   2^ 

pUm  D 


Cd  =-2B,  (19) 


substituting  from  equations  (9)  and  (13)  yields 


0     pUra  D 


^  ^^l  "^   ^"  (20) 


If  the  coefficients  do  vary  with  the  phase   a,   the  values  given  by 
equations  (18)  and  (20)  are  weighted  averages.   With  this  possibility 
in  mind  it  is  preferable  to  adopt 


■X—  =   Aj  sin  a  +  B   cos  a  |  cos  a  |  +  error  (21) 


or 


pUm  D 


F     TT  _    Da   .      Cd  I        1        ,  .„„v 

=  7"  Cm  •  —  sm  a r   cos  a   cos  a  +  error  (22) 


T^  ^  TT  2 

pUm  D  Um 


where  A  ,   B^  ,   Cm  and  Cd  are  constant  and  error  has  the  value 


error  =  A„  sin  3a  +  A  sin  5a 


+  B^   cos  3a  +  Bj.   cos  5a 


This  error  is  obtained  by  subtracting  the  computed  value  of  A-  sin  a 

I  n 

and  B    |  cos  a  |  cos  a  from  the  observed  F/pUm  D. 


13 


The  foregoing  formulation  was  developed  by  Keulegan  and  Carpenter 
[3]  and  equations  (18) ,  (20)  and  (22)  were  the  basis  for  the  computer 
data  reduction  program  developed  in  Appendix  A  and  presented  in 
Appendix  B. 


14 


III.   APPARATUS 

A.  TEST  BASIN 

A  test  basin,  16  feet  long,  4  inches  deep,  and  16.5  inches  wide  was 
constructed  from  3/4  inch  AC  marine  plywood.  The  side  walls  were 
supported  at  approximately  20  inch  intervals  with  1.5  inch  thick  plyvood 
gussets  and  two  2  by  6  inch  pine  stringers  along  the  top  edge  of  the 
tank.  Figure  1. 

The  presence  of  considerable  width-tolerance  variance  when  the  tank 
was  filled,  necessitated  the  addition  of  three  2  by  1  inch  channel 
braces,  welded  in  an  inverted  U-shape,  placed  over  the  top  of  the  channel, 
The  brace  design  allowed  for  unobstructed  motion  of  the  dynamometer 
along  the  length  of  the  tank,  an  advantage  not  afforded  by  structural 
cross  braces  along  the  tank  top.  Figure  2. 

All  joints  were  sealed  with  RTV  and  inside  surfaces  were  treated 
with  five  coats  of  penetrate  and  five  coats  of  epoxy  paint.   Two  3/4 
inch  thick  plexiglass  windows.  Figure  3,  (20  inches  by  12  inches)  were 
installed  near  the  tank  midsection. 

B.  DRIVING  MECHANISM 

The  driving  mechanism  consisted  of  a  "Dynamatic"  electric  motor  with 
an  eddy-current  clutch  and  Mark  III  controller.   This  system  utilized  a 
closed  loop  feedback  to  linearize  and  stabilize  its  shaft  speed.   The 
basic  control  system  used  is  presented  schematically  in  Figure  4. 
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Figure  2:      Cross  Braces 


■^Xjgf 


Figure  3:   Test  Section  Windows 
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Reference 
Input 

Error 
Detector 

•»■ 

Controller 

Controlled 

•-> 

Quantity 

Measure  of 

Controlled  C 

Quantity 

Figure  4:   Feedback  Loop 

More  specifically,  the  controlled  quantity  is  the  shaft  speed,  the 
reference  input  is  a  potentiometer  and  the  error  detector  is  a  summing 
amplifier.   The  controller  is  a  combination  of  a  power  amplifier, 
motor  and  eddy-current  clutch.   The  measure  of  shaft  velocity  is 
accomplished  by  means  of  a  tachometer  generator,  producing  an  output 
voltage  proportional  to  shaft  velocity. 

The  summing  amplifier  detects  differences  between  the  desired  speed 
and  actual  speed.   This  difference  is  fed  to  the  power  amplifier  which 
modulates  the  clutch  coil  current.   If  the  speed  is  too  low,  the  coil 
current  is  increased,  causing  the  eddy-current  clutch  to  transmit  more 
torque  to  the  load.  The  result  is  an  acceleration  of  the  shaft  to 
desired  speed. 

The  above  speed  regulation  operates  for  positive  loading  only. 
Since  the  minimum  coil  current  was  zero,  negative  loads  caused  free 
rotation  of  the  output  shaft. 

Both  positive  and  negative  loads  were  experienced  by  the  output 
shaft  due  to  the  connecting  rod's  weight  and  carriage  inertia.  This 
lead  to  uneven  output  shaft  speed,  particularly  at  low  RPM,  and 
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necessitated  the  application  of  a  continuous  positive  load.   This  load 
was  supplied  through  a  3  inch  wide  17  inch  diameter  brake  drum,  affixed 
to  the  face  plate,  with  a  120°  leather  lined  brake  shoe  riding  on  it. 

The  assembled  motor,  face  plate  and  brake  were  mounted  on  a  separate 
platform.  Figure  5,  to  isolate  any  inherent  mechanical  vibrations  from 
the  tank  structure.   A  7  foot  8  inch  long  by  1.5  inch  diameter  thin 
walled  aluminum  tube  was  then  used  to  connect  the  driving  mechanism  to 
the  carriage.   This  long  driving  rod  produced  a  near  sinusoidal  motion 
of  the  carriage. 

C .   CARRIAGE 

The  carriage  as  shown  in  Figure  6  provided  a  movable  mounting  plat- 
form which  was  positioned  over  the  test  tank. 

The  carriage  consisted  of  a  U-shaped,  welded  aluminum  structure 
fitted  with  four  3/4  inch  linear  bearings  located  at  the  ends  of  the 
hollow  cylindrical  side  members.   The  sides  of  the  platform  were 
constructed  of  two  inch  aluminum  pipe  fitted  with  aluminum  bearing  blocks 
at  each  end.   The  structural  transverse  member  of  the  carriage  was  a  2 
by  3  inch  aluminum  box  section  welded  to  the  end  of  each  pipe.   At  the 
opposite  ends  of  the  pipes  the  bearing  blocks  were  connected  by  a  3/8 
inch  diameter  transverse  shaft  to  accommodate  the  vertical  struts  which 
supported  the  cylinder. 

The  U-shaped  platform  was  mounted  through  linear  bearings  on  two 
case  hardened  3/4  inch  shafts.  The  shaft  ends  were  supported  by  two 
2  by  6  inch  steel  channels  that  transversed  the  width  of  the  tank.  These 
channels  were  clamped  to  the  stringers  along  the  top  of  the  side  walls. 
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Figure  5:   Driving  Mechanism 
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The  movable  part  of  the  carriage  consisted  of  two  long  arms  pivoted 
near  their  centers  and  mounted  on  a  drum  as  shown  in  Figure  7.   These 
arms  were  designed  to  hold  the  test  cylinders  at  a  maximum  centered  depth 
of  18  inches  submergence.   The  leading  edges  of  the  tapered  1/16  inch 
thick  polished  support  arms  were  sharpened  to  reduce  drag  and  experimental 
errors. 

A  counter  balance  was  placed  on  the  support  arms  opposite  the  test 
cylinder  and  adjusted  so  that  the  center  of  gravity  of  the  system  was  at 
the  pivot  point.   As  a  result  of  this  balance  adjustment  the  only  force 
(or  moment)  acting  on  the  arms  was  due  to  the  fluid;  all  other  inertial 
loads  were  balanced  out. 

A  third  arm,  half  the  length  and  perpendicular  to  the  supporting 
arms,  was  welded  to  the  center  drum.'  This  arm  served  as  a  2:1  mechanical 
amplifier  of  force  and  converted  the  torque  about  the  rotational  center 
of  the  drum  to  a  linear  force.   This  force  acted,  through  a  pinned  ball 
bearing  linkage  to  a  small  cantilever  beam  equipped  with  strain  gages. 
Figure  8. 

D.   CANTILEVER  BEAMS 

Two  cantilever  beams  were  constructed.  Figure  9; 

1.  0.25  inch  wide  by  0.35  inch  high  with  action  arms  of 
2.5  inch  and  4.5  inch 

2,  0.6  inch  wide  by  0.45  inch  high  with  action  arms 
1.5  inch  or  3.5  inch 

,Two  attachment  points  were  provided  on  each  beam  so  that  a  large  range 
of  forces  could  be  measured  without  overstressing  the  beam.   Stiffness, 
and/or  sensitivity,  could  be  adjusted  to  reduce  natural  frequency  vibra- 
tory motion  by  shortening  the  attachment  point  or  by  a  change  of  beams. 
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Four  strain  gauges  were  applied  to  each  beam  to  provide  a  means  of 
measuring  forces  imposed  upon  the  cylinder  and  transmitted  mechanically 
through  the  balance  assembly  to  the  beam.  .  •■ 

E.   TEST  CYLINDERS 

Three  cylinders  as  shown  in  Figure  10  were  constructed  to  provide 
maximum  flexibility  in  testing  over  a  wide  range  of  frequency  and 
displacement.   All  cylinders  were  fabricated  from  plexiglass  and  slotted 
to  accommodate  one  end  of  the  vertical  support  arms.   Bulkheads  were 
placed  adjacent  to  the  slots  to  provide  attachment  points  for  the  ends 
of  the  arms.   Silicone  sealent  was  used  to  prevent  leakage  after  the 
arms  were  attached.   This  allowed  the  cylinder  to  be  unflooded  resulting 
in  a  minimum  fixed  mass.   End  pieces  as  shown  in  Figure  11  were  also 
constructed  of  plexiglass.   The  ends  were  fitted  into  the  cylinder  with 
1/8  inch  "O"  rings  to  provide  the  necessary  seal.  The  outer  surfaces 
of  the  end  pieces  had  a  labyrinth  seal  cut  into  the  face  and  1/16  inch 
plexiglass  end  plates,  one  inch  larger  than  the  cylinder.   The  labyrinth 
seal  and  end  plates  were  provided  for  purposes  of  reducing  the  leakage 
around  the  ends  of  the  cylinder. 
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Figure  11;   Cylinder  End  Pieces 
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IV.   TEST  PROCEDURE 

A.  SYSTEM  SELECTION 

Primary  consideration  was  given  to  the  selection  of  the  cantilever 
beam  to  be  used  with  each  of  the  test  cylinders,  trade  offs  being  stiff- 
ness and  sensitivity.   It  was  concluded  that  the  4  inch  and  2.75  inch 
diameter  cylinders  would  be  run  on  both  beams  and  the  1.5  inch  diameter 
cylinder  tested  only  on  the  more  sensitive  of  the  two  beams.   This 
procedure  yielded  force  traces  that  could  be  cross  correlated  for  the 
4  inch  and  2.75  inch  cylinders  with  the  better  trace  utilized  for  data. 

Once  components  were  selected,  they  were  assembled,  balanced,  aligned 
to  the  tank  axis  and  clamped  to  the  tank.  edge.   The  system  was  then  finely 
balanced  in  the  dynamic  mode  by  adjusting  the  lever  arm  of  the  counter 
balance  weight  after  observing  strip  chart  traces  of  high  frequency 
oscillations  in  air.  The  counter  balance  weight  was  adjusted  so  that  no 
force  was  observed  when  oscillating  the  carriage  in  air.   The  tank  was 
then  filled  to  a  level    4  to  6  diameters  above  the  test  cylinder  and 
the  desired  amplitude  of  oscillation  and  frequency  set.  Amplitude 
settings,  i.e.,  connecting  rod  eccentricity  on  the  face  plate,  were  in 
half  diameter  increments.  The  frequency  setting  was  either  random  when 
observing  the  Reynolds  number  parameter  or  adjusted  to  maintain  a  given 
Reynolds  number  for  studies  of  the  amplitude/diameter  parameter. 

B.  SYSTEM  CALIBRATION 

1.   Force  Calibration 

A  length  of  fine  nicrome  wire  was  attached  to  the  test  cylinder, 
fare  lead  over  a  pulley  arrangement,  attached  to  the  top  of  the  tank. 
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After  a  half  hour  warm  up  period,  the  four  arm  strain  gage  bridge  was 
electrically  balanced  in  an  unloaded  state.   Then  accurate  weights  were 
hung  from  the  wire  and  the  resulting  deflection  of  the  recorder  pen  noted. 
The  strain  gauge  bridge  and  recorder  were  calibrated  in  both  the  positive 
and  negative  direction.   The  resultant  deflections  of  the  recorder  pen 
were  accurately  plotted  against  the  force  imposed.  Figure  12.  The  force 
calibration  proved  to  be  very  linear. 

This  calibration  was  conducted  prior  to  all  test  runs. 

2.      Event  Marker  Calibration 

Timing  of  the  moment  of  maximum  velocity  was  determined  by 
timing  contacts  on  the  drive  wheel.   The  ability  of  the  event  marker  to 
record  the  frequency  of  the  imposed  motion  and  the  phase  angle  between 
the  force  and  motion  was  dependent  upon  the  accuracy  of  the  chart  speed. 
However,  this  speed  was  checked  against  internal  and  external  timers  and 
no  discrepancy  was  found. 

C.   DATA  RUN 

The  system  as  selected  was  set  in  motion  and  recorder  started.   To 
avoid  wave  reflections  from  the  ends  of  the  tank,  the  system  was 
oscillated  for  less  than  ten  cycles  and  secured.  No  appreciable  standing 
wave  was  generated  under  this  procedure  and  the  system  was  allowed  to 
remain  at  rest  for  fifteen  minutes  before  a  new  data  set  was  initiated. 
After  each  data  run,  the  force  calibration  and  event  marker  position 
were  checked. 
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D.   STRIP  CHART  READ  OUT 

A  mean  curve  was  fared  through  the  force  recording  to  smooth  out 
natural  frequency  fluctuations  and  bearing  noise.  This  curve  was  then 
divided  into  twenty  divisions,  Figure  13,  and  mid-division  force  values 
were  read  and  correlated  to  its  respective  t/T  for  inputs  to  a  data 
reduction  computer  program,  Appendix  B . 
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Figure  13:   Sample  Force  Read  Out 
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V.   RESULTS 

The  purpose  of  the  experimental  program  was  to  study  the  effect  of 
the  parameter  A/D  on  the  added  mass  and  drag  coefficients  at  a  fixed 
Reynolds  number.   Also  a  limited  number  of  runs  were  conducted  to  deter- 
mine the  effect  of  Reynolds  number  on  these  coefficients  at  a  fixed  value 
of  relative  displacement,  A/D.  This  latter  series  of  runs  was  particularly 
significant  since  in  the  past  it  has  been  taken  for  granted  that  little 
or  no  Reynolds  number  effect  occurs  in  flows  of  the  present  type. 

Figures  14b  and  15b  show  the  total  calculated  force,  based  on 
Figures  14a  and  15a,  and  the  observed  force  for  two  typical  test  runs. 
These  figures  show  that  the  representation  of  the  total  force  as  the  sxim 
of  the  drag  and  inertia  components  is  adequate  for  practical  application. 
Figure  14  corresponds  to  a  case  dominated  by  inertia  while  Figure  15 
corresponds  to  a  larger  value  of  A/D  where  the  drag  component  is  the 
greatest.   The  error  computed  in  Appendix  A  is  the  difference  in  the 
observed  and  calculated  values  plotted  in  Figures  14  and  15. 

Figures  16  and  17  show  the  results  for  Cd  and  Cm  plotted  as  a 
function  of  A/D  for  a  smooth  and  rough  cylinder.   Also,  Keulegan  and 
Carpenter's  mean  value  line  [3]  is  plotted  for  comparison.   Keulegan 
and  Carpenter  conducted  their  studies  with  stationary  cylinders  and 
standing  waves  and,  therefore,  had  no  independent  means  of  controlling 
Reynolds  number  for  a  .fixed  A/D.   In  general,  however,  their  studies 
corresponded  to  fairly  low  Reynolds  number. 

It  is  noted  that  the  mean  Cm  line  taken  from  their  work  has  been 
reduced  in  value  by  1.0  to  account  for  differences  in  test  procedure 
(stationary  fluid  versus  moving  fluid) . 
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It  can  be  seen  that  Cd  values  which  correspond  to  rather  high  Reynolds 
number,  fall  below  Keulegan  and  Carpenter's  results.   However,  the  series 
runs  which  were  made  at  fixed  A/D  with  increasing  period  (decreasing 
Reynolds  number)  show  that  the  values  tend  toward  the  low  Re3molds  number 
results  of  Keulegan  and  Carpenter.   It  is  surprising,  however,  that  the 
effect  of  roughness  is  to  increase  the  drag  coefficient  so  as  to  approach 
Keulegan  and  Carpenter's  mean  line. 

Figure  17  shows  the  variation  of  Cm  with  A/D  along  with  the  results 
of  Keulegan  and  Carpenter.   It  is  noted  that  the  results  for  rough 
cylinders  and  for  higher  Reynolds  numbers  encountered  in  this  study  fall 
above  those  of  Keulegan  and  Carpenter.   It  is  interesting  that  roughness 
tends  to  have  little  effect  on  Cm  in  the  present  study. 

The  results  for  Cm  and  Cd  plotted  against  A/D  is  presented  in 
Figures  18  and  19,  respectively.   No  effort  was  made  in  these  plots  to 
keep  the  Reynolds  number  constant.   As  a  result  the  Cd  plot  shows  a 
large  amount  of  scatter;  the  Cm  plot,  however,  appears  to  correlate 
fairly  well  without  regard  for  the  Reynolds  number. 

Figure  20  and  21  show  the  added  mass  and  drag  coefficients  plotted 
against  A/D  for  a  fixed  Reynolds  number.   It  is  noted  that  the  scatter 
caused  by  the  non-constant  Reynolds  number  has  been  suppressed,  and  the 
results  correlate  well. 

Keulegan  and  Carpenter's  mean  lines  also  have  been  plotted  for 
comparison  on  Figures  20  and  21.   The  magnitudes  of  individual  points 
do  differ  but  the  same  trends  are  prevalent  throughout.  This  difference 
is  attributed  to  experimental  differences  (particularly  Reynolds  number) 
between  studies  and  inherent  experimental  error.   The  critical  region  of 
A/D  remains  in  the  region  of  2.25  to  2.5. 
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Figures  22  and  23  show  the  effect  of  Reynolds  number  on  the  added 
mass  and  drag  coefficients,  respectively.   Figure  21  shows  that  the  added 
mass  coefficient  is  almost  independent  of  Reynolds  number  and  highly 
dependent  upon  A/D  in  the  range  tested.   However,  Figure  23  shows  that 
Cd  is  strongly  affected  by  the  Reynolds  number  and  generally  it  shows  a 
decrease  with  increasing  Reynolds  number. 

Keulegan  and  Carpenter  [3]  state  that  apparently  no  correlation 
exists  between  the  two  coefficients  and  Re3molds  number.   However,  on 
plotting  their  data  in  the  format  of  Figures  22  and  23,  the  same  general 
trends  are  observed;  Cm  appears  to  be  independent  of  Reynolds  number 
while  Cd  is  strongly  dependent.   This  data  is  plotted  in  Figures  24  and 
25.   It  appears  that  Keulegan  and  Carpenter  overlooked  this  rather 
significant  feature. 
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VI.   CONCLUSIONS 


1.  Total  force  acting  on  a  cylinder  oscillating  harmonically 
is  comprised  of  two  components,  drag  and  inertia. 

2.  The  total  force  may  be  mathematically  modeled  using 
coefficients  of  added  mass  and  drag. 

3.  Cm  is  primarily  a  function  of  relative  displacement  A/D 
and  appears  to  be  nearly  independent  of  Reynolds  number. 

4.  Cd  is  a  function  of  both  relative  displacement  and  Reynolds 
nuniber . 


47 


APPENDIX  A 

A.   FORMULATION  OF  THE  COMPUTER  PROGRAM  FOR  DATA  REDUCTION 

A  computer  program  was  used  to  reduce  the  data.  The  purpose  of  the 
program  was  to  compute  Cm,  the  inertia  coefficient,  and  Cd,  the  drag 
coefficient.   The  starting  point  for  the  problem  was  Keulegan  and 
Carpenter's  [3]  formulation  for  these  parameters.   Specifically: 

2UmT   C    F  sin  a 


and 


Cm  =  -2 —  J  2 ^  "  (^) 

TT  D  0     pUm  D 


Cd  =  -3/4   S''  ^-^f^     da  (2A) 

0     pUm  D 


The  above  integrals  were  transformed  into  the  following  forms: 


Cm  =  -^ Y   F  sin  (^)  /(t/T) 

TT  D  L  p  A  ^-' 

2 

8-  p  DttLA  '-' 


where  A  denotes  the  amplitude  of  the  oscillatory  motion  and  (/"(t/T) 
was  conveniently  chosen  as  0.05. 

The  corresponding  error  for  each  data  set  was  computed  by  the 
following  relationship: 


(3A) 


(4A) 


A8 


error  =    [ ^^T"    7  F  sin  ^  /(t/T)]    sin    (^) 

2   p  DttVl    ^-^  ^  ^ 


2 

+    [ ^^-y-    y   F  COS    (^)   cTCt/T) 

2   p   Dtt  A  L 

-  0.02    y    F   cos    (~^)    cTCt/T)]    cos    (6tt   t/T) 


(5A) 


where : 

T  =  period  in  seconds 

D  =  cylinder  diameter  in  feet 

L  =  cylinder  diameter  in  feet 

t  =  time 
This  error  function  is  obtained  by  subtracting  the  computed  values  from 
the  observed. 

The  program  employed  time  increments  of  1/20  T.   where  t/T  takes  on 
values  from  zero  to  one.   For  each  run  at  constant  period  it  was  possible 
to  fix   t/T   increments  for  each  data  set. 
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B.   FLOW  CHART  OF  COMPUTER  PROGRAM 


DEFINE 
PARAMETERS 


1 

INITIALIZE 

1 

CONSTANTS 

y 

f 

INITIATE 
DATA 

SETS 

\ 

■ 

READ 

PARAMETERS 

.  ,> 

/ 

"s 

COMPUTE 
COEFFICIENTS 

z,  z^  z. 

12   3 

> 

' 

COMMENCE 

DATA 

REDUCTION 

« 

f 

READ 

DATA 

SETS 

^ 

f 

COMPUTE 

TRIG 

ARGUMENT 

t 

f 
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i  '          -^ 

i 

COMPUTE 
FORCE 

\r 

PRINT 
COMPONENT 

^  f 

COMPUTE 
Cd  Cm 

U 

COMPUTE 
ERROR 

^ 

' 

PRINT 
ERROR 
Cd  Cm 

5 

'' 

■ 

STOP 
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APPENDIX  B 


COMPUTER  PROGRAM 

10    FORMAT(3F10.4,I10) 

15    FORMAT {  •  1' , lOX,  •DIA=«, F8.4,5X, • AMP=« ,F8.4,5X,«PER  =  ',F8 
1.4,5X»'Z1=' 

1F8.4»5X» •Z2=', F8.4»5X,'Z3=« ,F8.4,5X,'NCARD=',I3) 
1«»8X,«F' ,11X 
20    FORMAT (• J' , 3X ♦ ' T IM E/PER* ,7X,' ALPHA* ,7X, 'C0SA',8X, "SINA 

1, 'FCOSA '.SX, 'FSINA' ) 
25     FORMAT  (2F10. 4) 
30    FORMAT( • 0*  ,7FL2.4} 

35    FORMAT ( '0»  WXt  •CM=« ,7X,  '00=' »7X, •BETA=« ) 
40    FORMATCO*  ,3F12.4) 
45    FORMAT (  ' 0*  ,  lOX , » A3= • »F 8.4 , 5X » • Bl  =  '  , F8 .4 » 5X  ,  • B3  =  • , F8 .4 , 

15X, • B3P=', F8.4) 
50    FORMAT (• 0'  t3X ,' T I  ME/PER'  ,7X,«  ALPHA' »7X,  'S IN3A»,7X,  'COS 

13A«,7X, 'ERROR' ) 
55     FORMAT  CO'  »5F12.4J 
C  TIME=DIMENSIONLESS    TIME     (TI ME/ PERIOD) 

C  BETA=DIMENSIONLESS    DISPLACEMENT    {  UMAXs^PER/DI  A) 

C  L=CYLINDER    LENGTH    IN    FEET 

C  DIA=DIAMETER    OF    CYLINDER    IN    FEET 

C  AMP=AMPLITUDE    CF    MOTION    IN    FEET 

C  PER=PERIOD    IN    SECS     (CHART    PERIOD/CHART    SPEED) 

C  ZI=FORCE    COEFFICIENTS 

C  CM=INERTIA    COEFFICIENT 

C  CD=DRAG    COEFFICIENT 

C  ERROR=REMAINDER    FUNCTION 

C  A1,BI=F0URIER    COEFFICIENTS 

C  K=CONSTANT 

C  N=NUMBER    OF    DATA    SETS 

G=32.174 
RH0=62.4/G 
K=0.2 

PI=3. 14159 
L=1.375 
C  INITIATE    DATA    SETS 

N=9 

DO  100  1=1, N 
A3=0.0 
B1=0.0 
B3=0.0 
CM=0.0 
CD=0.0 
C  READ    IN    PARAMETERS    WHICH    CHARACTERIZE    THE    DATA    SET 

READ{5tl0)    DIA,AMP,PER,NCARD 
C  COMPUTE    FORCE    COEFFICIENTS 

Zl=2*PER**2/(PI«=^3-0IA**2'i=L*RH0*AMP) 
Z2=-3^PER**2/{  8*RH0-DI  A'!'PI=^L^AMP**2) 
Z3=PER**2/(2*RH0^DI  A^PI*=*2=^  AMP**2*L  ) 
C  COMPUTE    BETA 

BETA=2*PI*AMP/DIA 
C  PRINT    OUT    PARAMETERS    AND    COEFFICIENTS 

WRITE(6,15)     DIA,AMPtPER,ZltZ2,Z3,NCARD 
C  PRINT    OUT    FORMAT 

WRITE(6,20) 
C  COMMENCE    DATA    REDUCTION 

DO    1000    J=1,NCARD 
READ(5,25)    TIMEtF 
ALPHA=2=<=PI*TIME 
SINA=SIN(ALPHA) 
SIN3A  =  SIN(3'^  ALPHA) 
COSA=COS(ALPHA) 
C0S3A=C0S(  3-ALPHA) 
FSINA=0.05*F*SINA 
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FSIN3A=0.05=^'F*SIN3A 

FCOSA=0.05-F=^COSA 

FC0S3A=0.0  5=^F*C0S3A 
C  SUM    COMPONENTS 

CM=FSINA+CM 

CD=FCOSA+CD 

A3=FSIN3A+A3 

B1=FC0SA+B1 

B3=FC0S3A+B3 
C  PRINT    OUT    COMPONENTS    FOR    DRAG    AND    INERTIA    COEFFICIENT 

WRITE(6,30)    TIME, ALPHA, COSA,SINA,F,FCOSA,FSINA 
1000    CONTINUE 
C  COMPUTE    COEFFICIENTS 

CM=Z1*CM 

CD=Z2*CD 

A3=Z3*A3 

B1=Z3*B1 

B3P  =  Z3'!'B3-K*B1 
C      PRINT  OUT  INERTIA 

WRITE{6,35) 

WRITE(6,40)     CM, CD, BETA 
C  COMPUTE    ERROR 

C  PRINT    OUT    ERROR    FORMAT 

WRITE(6,45)     A3,B1,B3,B3P 

WRITE(6,50) 

TIME=0.025 

DO    1500    K=1,NCARD 

ALPHA=2*PI*TIME 

SIN3A=SIN(3-ALPHA) 

C0S3A=C0S( 3=ALPHA) 

ERR0R  =  A3*S  IN3A  +  B3P=^C0S3A 
C  PRINT    OUT    ERROR    COMPONENTS    AND    ERROR 

WRITE (6, 55)    TIME,ALPHA,SIN3A,C0S3A,ERR0R 

TIME=TIME+0.05 
1530    CONTINUE 
100    CONTINUE 

STOP 

END 


AND  DRAG  COEFFICIENTS 
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APPENDIX  C 

A.      CALCULATION  OF   INERTIAL  AND   DRAG  FORCES 
Given 

T,   Cd,   Cm,   D,   Am,      and 

p   =  1.93   slugs 


F  =  Fj  +  Fjj  =  Cm   p  Ao  ^     +  j  Cd  D   p  U    |u 


Therefore; 


Therefore; 


r^      n    A      dlj  _      ottD^     dU 

F^.   =   Cm   P   Ao  -r-     =      Cm    ^    .        -j- 

I  dt  4        dt 


U  =   -Um  cos   o   t         -jz  =  Um  a   sin  a   t 

at 


a     =  2  tt/T 
Um  =  2   TT  — 


dU        2tt  Am     27f 


Att     Am        , 

sin  a   t 


i^ 


2  2 

_          ^     pTT   D        Att     Am        ,      2TTt 
Fj  =  Cm  — ^    —     sm  ^- 


5A 


Let  a  =  ^  =  a(t) 


with  a  =  27T     t/T 


3     2 

^  Cm  P    TT     D     Am 

r,  =  X sm  a 

T 


Fjj  =  Y  Cd  D  P    U    |u| 


substituting   for  U 


1  2    I  I 

F       =  -  y  Cd  D  p  Urn      |    cos  a   t    |    cos   a   t 


substituting  for  Um 


therefore 


2   2 

1  _, ,  ^   Att  Am  I 

F   =  -  p-  Cd  D  p —   cos  a   t    \    cos  a  t 


2        2 

„      „,   2D  p  TT  Am  I 

F_  =  -  Cd cos   a   I  cos   a 


when  a  is  substituted  for  at. 
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CALCULATED  DRAG  AND  INERTIA  FORCE  FOR: 


T  =  1.60  seconds 


Cm  =  1.1722 


Cd  =  0.7947 


p  =  1.93  slugs 


A  =  0.1667  ft. 


D  =  0.3333  ft, 


F  =  0.508  sin  a 

F  =  -0.110  cos  a  I  cos  a 


F  =  F  + 

c  -^i 

^D 

Time 

Inertial 

Drag 

Calculated 

Observed 

Period 

Force 

Force 

Force 

Force 

0.25 

0.079 

-0.107 

-0.027 

0.05 

0.75 

0.230 

-0.087 

0.143 

0.20 

1.25 

0.360 

-0.055 

0.304 

0.33 

1.75 

0.453 

-0.023 

0.430 

0.43 

2.25 

0.502 

-0.003 

0.499 

0.50 

2.75 

0.507 

0.003 

0.505 

0.50 

3.25 

0.453 

0.023 

0.475 

0.50 

3.75 

0.360 

0.550 

0.414 

0.45 

4.25 

0.230 

0.087 

0.318 

0.35 

4.75 

0.080 

0.107 

0.187 

0.17 

5.25 

-0.080 

0.107 

0.048 

0.05 

5.75 

-0.230 

0.087 

-0.143 

-0.10 

6.25 

-0.360 

0.055 

-0.304 

-0.25 

6.75 

-0.453 

0.023 

-0.430 

-0.40 

7.25 

-0.502 

0.003 

-0.499 

-0.45 

7.75 

-0.502 

-0.003 

-0.505 

-0.50 

8.25 

-0.453 

-0.023 

-0.475 

-0.50 

8.75 

-0.360 

-0.055 

-0.414 

-0.45 

9.25 

-0.230 

-0.087 

-0.318 

-0.35 

9.75 

-0.079 

-0.107 

-0.187 

-0.10 
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CALCULATED  DRAG  AND  INERTIA  FORCES  FOR: 

T  =  1.39  seconds      Cm  =  0.9838         Cd  =  1.8856 

p  =  1.93  slugs        ^  ==  0.3437  ft.      D  =  0.2991  ft, 


Fj.  =  0.667 

sin  a 

Fj,  =  -3.17 

cos  a   ( 

cos  a 

X 

F^  =  Fj  +  1 

^D 

Time 

Inertial 

Drag 

Calculated 

Observed 

Period 

Force 

Force 

Force 

Force 

0.025 

0.087 

-0.983 

-0.896 

-0.900 

0.075 

0.252 

-0.800 

-0.548 

-0.600 

0.125 

0.393 

-0.504 

-0.111 

-0.300 

0.175 

0.496 

-0.207 

0.289 

0.000 

0.225 

0.550 

-0.024 

0.525 

0.300 

0.275 

0.550 

0.024 

0.574 

0.600 

0.325 

0.496 

0.207 

0.703 

0.900 

0.375 

0.393 

0.504 

0.897 

1.200 

0.425 

0.252 

0.800 

1.052 

0.900 

0.475 

0.087 

0.983 

1.070 

0.700 

0.525 

-0.087 

0.983 

0.896 

0.900 

0.575 

-0.252 

0.800 

0.548 

0.700 

0.625 

-0.393 

0.504 

0.111 

0.250 

0.675 

-0.496 

0.207 

-0.289 

-0.250 

0.725 

-0.550 

0.024 

-0.525 

-0.600 

0.775 

-0.550 

-0.024 

-0.574 

-0.800 

0.825 

-0.496 

-0,207 

-0.703 

-1.100 

0.875 

-0.393 

-0.504 

-0.897 

-1.100 

0.925 

-0.252 

-0.800 

-1.052 

-0.700 

09.75 

-0.087 

-0.983 

-1.070 

-0.650 
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APPENDIX  D 


DATA  RUNS  FOR  REYNOLDS  NUMBER  2.7  x  IQ 


Runs 


Cyl.  Dia. 


A/D 


Cm 


Cd 


Rl 

1.5  in. 

1.5 

0.46 

1.4692 

1.7197 

R2 

1.5 

2.0 

0.60 

0.2134 

2.1625 

R3 

1.5 

2.5 

0.78 

-0.1047 

2.4826 

R4 

1.5 

3.0 

0.91 

0.3988 

2.0726 

R5 

1.5 

3.5 

1.05 

0.2462 

2.1140 

R6 

1.5 

4.0 

1.21 

0.7225 

2.0300 

R7 

1.5 

4.5 

1.37 

0.0249 

2.0825 

R8 

1.5 

5.0 

1.56 

-0.3539 

2.2148 

R9 

2.75 

0.75 

0.75 

1.4791 

1.1915 

RIO 

2.75 

0.75 

0.75 

1.4594 

1.0918 

Rll 

2.75 

1.0 

1.0 

1.2939 

1.1664 

R12 

2.75 

1.0 

0.96 

1.2832 

0.9221 

R13 

2.75 

1.5 

1.50 

0.8159 

1.7550 

R14 

2.75 

1.5 

1.45 

0.8762 

1.9681 

R15 

2.75 

2.0 

2.06 

-0.6060 

2.9893 

R16 

2.75 

2.5 

2.53 

0.1127 

2.3252 

R17 

2.75 

2.5 

2.57 

-0.0085 

2.4444 

R18 

4.0 

1/4 

0.53 

1.4959 

1.6421 

R19 

4.0 

3/4 

1.50 

1.5350 

0.9485 

R2G 

4.0 

1/2 

1.06 

1.4632 

1.1187 

R21 

4.0 

1.0 

2.15 

1.3128 

1.1295 
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DATA  RUNS   FOR  VARIABLE   REYNOLDS   NUMBER 


Run 


A/D 


1. 

1/2 

0.84 

2. 

1/2 

0.60 

3. 

1/2 

0.48 

4. 

1/2 

0.40 

5. 

1/2 

0.33 

6. 

1 

0.83 

7. 

•  1 

0.61 

8. 

1 

0.48 

9. 

1 

0.39 

10. 

1 

0.33 

11. 

1.5 

1.39 

12. 

1.5 

0.83 

13. 

1.5 

0.62 

14. 

1.5 

0.48 

15. 

2 

1.30 

16. 

2 

0.85 

17. 

2.5 

1.32 

18. 

2.5 

0.85 

19. 

2.5 

0.55 

20. 

1/2 

1.26 

21. 

1/2 

0.83 

22. 

1 

1.36 

23. 

1 

0.86 

24. 

1.5 

1.35 

25. 

1.5 

1.05 

26. 

1.5 

0.84 

Runs   1  - 

-  19 

Runs  20  - 

-  26 

Cm 

1.4201 
1.4358 
1.4395 
1.4244 
1.2812 
1.2561 
1.3029 
1.3375 
1.3154 
1.2295 
0.9838 
0.9190 
1.0179 
0.9710 
0.8083 
0.7454 
0.5224 
0.7891 
0.9003 
1.4857 
1.2976 
1.4547 
1.3394 
1.1263 
1.0685 
1.1361 


2.75  inch  diameter  cylinder 
4.00  inch  diameter  cylinder 


Cd 

Re  X  10 

1.0556   ' 

1.61 

1.1032 

2.25 

0.7405 

2.82 

0.3118 

3.38 

1.2897 

4.10 

0.9318 

3.26 

0.7524 

4.43 

0.5872 

5.63 

0.7546 

6.39 

0.9274 

8.19 

1.8856 

2.92 

1.1842 

4.89 

1.0826 

6.54 

1.0334 

8.45 

1.3861 

4.16 

1.2833 

6.36 

1.6231 

5.12 

1.1180 

7.95 

0.8348 

12.30 

1.2957 

2.26 

1.6856 

3.43 

0.8183 

4.20 

0.8546 

6.64 

1.1111 

6.35 

1.0063 

8.17 

0.8950 

10.10 

-4 


59 


BIBLIOGRAPHY 


1.  Reid,  R.  0.  and  Bretschneider ,  C.  L.,  "The  Design  Wave  in  Deep  and 

Shallow  Water  and  Forces  on  Pilings  and  Large  Submerged  Objects", 
A.  and  M.  College  of  Texas,  Dept.  of  Oceanography,  Tech  Report  on 
Contract  NOy  27474,  DA  -  49  -  005  -  eng  -  18,  and  N7  onr  -  48704, 
36  pp.,  February  1954  (unpublished). 

2.  Weigel,  R.  L. ,  "Oceanographical  Engineering",  Prentice  -  Hall,  1964. 

3.  Keulegan,  G.  H. ,  and  Carpenter,  L.  H. ,  "Forces  on  Cylinders  and  Plates 

in  an  Oscillating  Fluid",  Journal  of  Research  of  the  National 
Bureau  of  Standards,  v.  60,  no.  5,   p.  423-440,  May  1958. 

4.  Heinzer,  A.  J.,  "Experimental  Investigation  of  the  Wake  of  an 

Oscillating  Cylinder",  Unpublished  Master's  Thesis,  University  of 
Houston,  Houston,  Texas,  1968. 

5.  Sarpkaya,  T. ,  and  Garrison,  C.  J.,  "Vortex  Formation  and  Resistance 

in  Unsteady  Flow",  Journal  of  Applied  Mechanics  Transactions  of 
the  ASME,  v.  30,  series  E,  no.  1,   p.  16-24,  March  1963. 

6.  Morison,  J.  R. ,  O'Brien,  M.  P.,  Johnson,  J.  W.,  and  Schaaf,  S.  A., 

"The  Force  Exerted  by  Surface  Waves  on  Piles",   Journal  of  Petroleum 
Technology,  American  Institute  of  Mining  and  Metallurgical  Engineers, 
V.  189,  p.  149-154,  May  1950. 

7.  McNown,  J.  S.,  and  Wolf,  L.  W. ,  "Resistance  to  Unsteady  Flow:   I. 

Analysis  of  Tests  with  Flat  Plate",  Engineering  Research  Institute 
of  Michigan,  2446  -  I  -  P,  June  1956.   Internal  report  to  Sandia 
Corporation. 


60 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Asst.  Professor  C.  J.  Garrison  3 
Code  59Gm 

Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Professor  T.  Sarpkaya  1 
Code  59Si 

Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Lt.  John  R.  Driscoll,  Jr.  USN  2 
877  Crestline  Drive 

Blue  Bell,  Pennsylvania  19422  • 

6.  Department  of  Mechanical  Engineering  1 
Library,  Code  59 

Naval  Postgraduate  School 
Monterey,  California  93940 


61 


iREPORT     TITLE 


I,   DESCRIPTIVE   NOTES  (-Type  ot  fporfndjnclui  iv»  dmtei) 

Master's  Thesis;  December  1972 

r.  AUTMOR(S)  (First  num»,  middlt  initial,  /*»«  nam») 


|.    REPORT    DATE 

December  1972 


UNCLASSIFIED 


Socvintv  Cla?.sific3tion 


DOCUMENT  CONTROL  DATA  -R&D 

^S.cun.y  c,.,..,,c.„on  ol  „„..    tody  o,  et„r.c,  .nd  /ndeWn,  .nno„„on  n.usr  6.  .nt.r.d  >vhen   ,He  o..r.„  report  /,  cU,.,n,d) 

(ifcuniy    ti«  j_^     ^^ ^    SECURITY    CLASSIFICATION 

ORIGINATING    ACTIVITY   (Corporml*  author) 


Naval  Postgraduate  School 
Monterey,  California  93940 


Unclassified 


26.  CROUP 


Forces  on  Cylinders  Oscillating  in  Water 


John  R.  Driscoll,  Jr. 


CONTRACT    OR    GRANT    NO. 


b.    PROJECT   NO. 


7«.     TOTAL    NO.    OF    PACES 

63 


7b.    NO.    OF    REFS 

7 


»a.    ORIGINATOR'S   REPORT   NUMBERHI 


9b.  OTHER   REPORT  NO(S)  (Ar^y  otharnumbara  that  may  ba  aaalgnad 
Ihit  raport) 


10      DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


II.    SUPPLEMENTARY    NOTES 


12.    SPONSORING   MILITARY    ACTIVITY 

Naval  Postgraduate  School 
Monterey,  California   93940 


13.  ABSTR  AC  T 


A  circular  cylinder  immersed  in  an  infinite  fluid  and  oscillated 
perpendicular  to  its  axis  is  acted  upon  by  a  fluid  dynamic  force.   Thxs  force 
is  generally  considered  to  be  composed  of  two  components,  a  drag  and  an  added 
mass  component.   The  relative  importance  of  these  are  functions  of  the  frequency 
and  displacement  of  the  cylinder. 

It  was  the  purpose  of  this  study  to  determine  experimentally  these  two 
components  as  a  function  of  the  ratio  of  the  displacement  amplitude  to  the 
cylinder  diameter,  and  the  Reynolds  number  based  on  the  maximum  velocity. 


DD;NrJ473 

S/N    0101 -807-681 1 


(PAGE     1  ) 


62 


UNCLASSIFIED 


Security  Cl««iific«tion 


l-ai408 


UNCLASSIFIED 


Security  Classifiration 


K  e V   woMot 


^ 


UNSTEADY  FLOW 

DRAG 

ADDED  MASS 


DD  ,r:..1473  'BACK) 


5/N    0101 -807-6S21 


ROLE  W  T 


ROLE  W  T 


UNCLASSIFIED 


63 


Security  CU»*ific«tion  a-31409 


Thes » s 


0753     -^^^lIes.oncv^^ 


osci 


Hil46G 


Thesis 

0753         Driscoll 
c.l  Forces  on   cylinders 

osci  1  lating    in  water. 


thesD753 

Forces  on  cylinders  oscillating  in  water 


3  2768  002  00690  0 

DUDLEY  KNOX  LIBRARY 


